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Abstract 

A stochastic linear transport equation with multiplicative noise is considered 
and the question of no-blow-up is investigated. The drift is assumed only 
integrable to a certain power. Opposite to the deterministic case where 
smooth initial conditions may develop discontinuities, we prove that a certain 
Sobolev degree of regularity is maintained, which implies Holder continuity of 
solutions. The proof is based on a careful analysis of the associated stochastic 
flow of characteristics. 
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1. Introduction 

Consider the stochastic linear transport equation in Stratonovich form 

du , „ „ dW n 

— + b ■ Vu + aVu o — - = 0, u\ t=0 = u . 

at dt 

Here W = {W t ) t>0 is a <i-dimensional Brownian motion defined on a filtered 
probability space J?, J^, P), the drift b : [0,T] x R d R d is a given 
deterministic vector field, uGl and uq : lR d — > R are given and the solution 
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u = u(x, t) will be a scalar random field on (O, J^, P) defined for (t, x) G 
[0, T] x R d . 

We deal with the problem of singularities of u starting from a regular 
initial condition uq. When o = and b is not Lipschitz, singularities may 
appear, in the form of discontinuities (or blow-up of derivatives), as in the 
simple example d — 1, b (x) = —sign(x) \/Jx\: any non-symmetric smooth 
initial condition w develops a discontinuity at x = for any t > 0, be- 
cause there are different, symmetric, initial conditions xq for the associated 
equation of characteristics 

x' it) = b(t,x(t)) , x(0) = x 

which coalesce at x = at any arbitrary positive time. Opposite to the 
question of uniqueness of weak L°° solutions, where positive results have been 



given under relatively weak assumptions on b, see for instance and it 
seems that good results of no blow-up are not available in the deterministic 
case when b is not Lipschitz. 

The purpose of this paper is to show that, for o ^ and b of class 

beL q p :=L q (0, T; L p (R d , R rf )) , 

d 2 

p,q>2, - + -<1 1 
V Q 

some regularity of the initial condition is maintained, in particular disconti- 
nuities do not appear. We prove the following result. 

Theorem 1. If cr ^ 0, (|7]j holds and uq G fl r >iH /1 ' r (M. d ) then there exists a 
unique solution u such that 

P{u(t,-)en r > 1 wt; r c {R d )) = l 

for every t G [0, T] . 

The unique solution in this class is given by a representation formula, in 
terms of Uo, involving a weakly differentiate stochastic flow. By Sobolev 
embedding theorem, u(t, •) is a- Holder continuous for every a G (0, 1), with 
probability one. Hence, from smooth initial conditions, discontinuities cannot 
arise. 

The precise formulation of the concept of solution and other details are 
given in the sequel. 
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The intuitive idea is that, opposite to the deterministic case, when cr^O 
the characteristics cannot meet. They satisfy the stochastic equation 



dX t = b (X t , t) dt + adW t (2) 

which generates, under assumptions ([I]), a stochastic flow of Holder continu- 
ous homeomorphisms, with some weak form of differentiability. The existence 



of an Holder continuous stochastic flow has been proved in [12|, [13|, |21l. A 
differentiability property in terms of finite increments has been given in [l3j . 
Here we establish Sobolev type differentiability. A similar Sobolev regularity 



of the flow is investigated in 20| by different tools (Malliavin calculus). See 
also [l9|. 

The assumption has been introduced in the framework of stochastic 



differential equations by [17| who have proved strong uniqueness. In the fluid 
dynamic literature, with < in place of <, it is known as the Ladyzhenskaya- 
Prodi-Serrin condition. One of its main consequences is that it gives uniform 
bounds on gradients of solutions to an auxiliary parabolic problem (see The- 
orem [2] below) essential for our approach, along with good properties of the 
second derivatives. 

The possibility that noise may prevent the emergence of singularities is 
an intriguing phenomenon that is under investigation for several systems. 
For linear transport equations with b 6 L°° ([0,T] ; (M d ,lR a! )) it may be 
deduced from [15| (the result presented here is more general). For nonlinear 
systems there are negative results,_like the fact that noise does not prevent 
shocks in Burgers equation, see [14| . and positive results for special kind 
of singularities (collapse of measure valued solutions) for the vorticity field 
of 2D Euler equations, see 16j], and for ID Vlasov-Poisson equation, see 



lo| . Moreover, for Schrodinger equations, there are several theoretical and 
numerical results of great interest, see [3[ - Q. We do not list here the results 
concerning the restored uniqueness due to noise and address to the lecture 
note [14J on this subject. 

After the result of Theorem [TJ it remains open the question whether the 
solution is Lipschitz continuous (or more) when uq G PI/ 1 ' 00 (lR d ) (or more). 
In dimension d we think that this is a difficult question under assumption 
The answer is positive when b e L°° ([0,T] ; C% (R d ,R d )) because the 
stochastic flow is made of diffeomorphisms, see [15j and it is also positive 
in dimension d = 1 for certain discontinuous drift b, including for instance 
b (x) = sign (x), see 
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It must be emphasized that, although this "regularization by noise" may 
look related to the regularization produced by the addition of a Laplacian 
to the equation, in fact it preserves the hyperbolic structure of the equation. 
The equations remain reversible and the solution at time t is, in the problem 
treated in this paper, just given by the initial condition composed with a 
flow. If the initial condition has a discontinuity, the solution also has a 
discontinuity; no smoothing effect is introduced. However, the emergence of 
singularities (shocks in our case) is prevented. 

The work is organized as follows. In Section [2] we present some results 
on regularity and approximation properties of the flow associated to the 
SDE (j5J). They are obtained via the study of an associated SDE and the 
regularity of its solutions. The main results are contained in Lemmas [3] and 
|5l while more technical results are collected in the Appendix (Section EJ). In 
Section [3] we define weakly differentiable solutions of the SPDE and prove 
their existence in Theorem [TOJ A technical result on convergence of random 
fields in Sobolev spaces is left to the last Appendix. Finally, uniqueness of 
weakly differentiable solutions of the SPDE is proved in Section HJ 

2. Convergence Results 

In this section we present some technical results on an associated SDE 
that we will study as an intermediate step to obtain some regularity and 
approximation properties of the flow associated to the SDE (J2J). The main 
results are contained in Lemmas [3] and [51 

Let us start by setting the notation used and recalling some results. We 
will use the following auxiliary SDE, introduced in (l3| : 



dY t = \U(t,%\x))dt + (VU{t,ju l {x)) + Id)dW t , Y = x . (3) 



The link between this SDE and the one presented in the introduction is given 
by the C 1 -diffeomorphism j t : Y t = j t o X t o where jt( x ) — x + U(t,x). 
Here, U : R d+1 -> R d is the solution of the PDE 




Hl p (T) := L q (0,T; W 2 ' p (R d )) D W 1 * (0,T; L p (R d )) ; 
we report here the precise result, given by [13| , Theorem 3.3]. 
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Theorem 2. Take p, q such that (QJ) holds, A > and two vector fields 
b,f(t,x) : IR d+1 — > M. d belonging to L^(T). Then in iff p (T) there exists a 
unique solution of the backward parabolic system 

d t u + |Aw + b ■ Vu - Am + / = ( . 

For this solution there exists a finite constant N depending only on d, p, q, T, A 
and such that 

h| p (t) < N\\f\\ Li{T) . (5) 



u 



We will use the result of this theorem with / = b. 

Let b n be a sequence of smooth vector fields converging to b in L q p . Let U 
be the unique solution to the PDE (J3J) provided by the above Theorem and 
U n the solutions obtained using the approximating vector fields b n . Lemma 
T2l shows that the vector fields U n converge in H\ v to U. 

In 13] is also proved the existence of Holder flows of homeomorphisms 
for the two SDEs above, which we denote by (f>t(-) for the SDE ([2]), and ipt(') 
for (E]). We will use to denote the flows obtained for the approximating 
vector fields b n , and for the flows corresponding to the auxiliaries SDEs 
obtained via the diffeomorphisms 7™ = Id + U n (t, ■). We will use (j)Q n (-), and 



ipQ (•) for the inverse flows. 

We can now state and prove the two main regularity results on the flows </>o' n . 



Lemma 3. For every R > 0, p > 1 and x,y G B 



lim sup sup E 

n ^°°te[o,T}xeB R 



^(x)-0^)| p <C p , T \x-y\ p 



In particular, 



lim sup sup E [|0o' n (x) — 0q (x)| p 1 = 0. (6) 

n ^°°te[o,T] xeB R 



Proof. Step 1 (preliminary estimates). By lemma [T2l we have that for every 
r > there exist a function / s.t. lim^oo f{n) = and 

1 



sup sup VC/(i, x)\ < 



2' 



xm d te[o,r] 
sup sup \U n {t,x) - U(t, x)\ < f{n) 

xeB r tG[0,T] 

sup sup \VU n (t,x) - VU(t,x)\ < f{n) 

x&B r te[0,T] 
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Since 4>t(x) is jointly continuous in space and time, there exist an r < oo 
s.t. the image of Br x [0, T] will be contained in B r for all t < T. In the 
following we will always take x,y £ Br. It follows that 

\U n (t,ff(x))-U(t,<j> t {y))\ <f(n) + ±\<%(x)-Mv)\ , 

\vu n (t,w(x))-vu(t,My))\ < f(n) + \vu n (t,w(x)) - vu n (t,<p t (y))\ ■ 

To shorten notation, we will write (fi n and (ft to denote <fit(x) and <fit(y), 
U n ((j) n ) and U n (x) to denote U n (t,(f) n ) e f/ n (0,x), etc. The same holds for 
the flows of the SDE ([3]). From the definition ip™ = j t o 0™ o (70 ) _1 and the 
properties of the diffeomorphisms 7™ obtained from Lemma [12] and Remark 
[T5| we immediately have 

iv n - > \r - + tn<n - u n m - /(n) 

>^I0"-0|-/H (7) 
2(V-Vl + /W) > I0 n -0| 

and 

\V-^\<\\r-^\ + f{n) ■ (8) 



Step 2 (computations). We start by proving the convergence of the flows 
of the auxiliary SDE ([3]). By Ito formula, for any a > 2 



id 

a 



0-2 



+ <(V n - V), (VfT(<T) - W(0)) • dWi), 
a - 1 



Tr [W 



vu 



0-2 



Let us analyze the three terms Ai, A 2 , A 3 . Using (J7J) we have 



A 1 <\\r-^\{f{n) + \\4> n -4>\) 



dt 



< A|V™ - Vfdi + 2Xf(n)\^ n - V| dt . 



Idt 
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Since VU n is bounded (uniformly in n, see Lemma IT!?]) and by (|26p \ijj n \ a 
belongs to L 2 (Q x [0,T]) for any a > 1, we can write A 2 = dM™, where for 
every n, dM™ is the differential of a zero mean martingale. As for the third 
term, using twice the inequality (a + (3) 2 < 2(a 2 + (3 2 ) and the estimates of 
the first step, we get 



d 2 (a - 1) 



A* < \ 2\VW 



< 



Vf/ n (0)| 2 + 2/ 2 (n)}dt 



*dA? + 2f 2 {n) dt 



< 8\ip n - ^| <L4£ + 8f 2 (n) dA™ + 2f 2 (n) dt 



where for every n 



A n t :=2 



\VU r ' 



VU r ' 



-l{4>n^4> s }ds 



(9) 



is a nondecreasing adapted stochastic process, with Aq = 0, and uniformly 
in n E[A%] < C < oo, see Lemma H Set := [Ad 2 a(a - 1)]A%. From 
the above estimates and after renaming M t (which remains a zero mean 
martingale), we get 



d(e- B "|r-Vr) <e" 



a\\ip n - Vl° + 2a\f(n)\ip n - ip\ 



a-l 



dt 



+ dM t + f 2 {n)e~ B ?\ip n - ip\ u *d£ t n 



ia-2 



+ d 2 a(a - 1) e~^f 2 (n)\ij n - ^| a ~ 2 dt . 

Integrating in time, taking the expected value, and finally the supremum 
over t G [0, T], we get 



sup E 

t€[0,T] 



e- B ?U?-ip t \ 



< < -Vo \ +a\E 



e~ B °\i) n s -i) s \ a ds 



+ C aAX f{n)E 



+ f(n)E 



BV: I |„/.n 



l a-l 



V£-iM + /(«) V£-ik 



ia-2 



e -B?U«_^ 



ia-2 



d£? 



ds 
(10) 



The expected value in the second line above is bounded uniformly in n. This 
fact is easily seen using for each term Holder inequality together with the 
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integrability properties of the flows ijj n and of the exponential of the processes 
B™, provided by fTStj]) and Lemma dH respectively. We claim that also the 
expected value of the last line is bounded. 



Claim 4. There exists a constant C s.t. for every n and p > 



E 



-B"\i.n 



o 



< C . 



Proof of the Claim. Using the definition of B™ we can rewrite the term 
on the left hand side as 



E 



-B 1 } L i.n 



. p \VU{s,<j>«)-VU{s, 



-ds 



'-jo 



Using Holder inequality, for some e > 1 small (to be fixed later) and k the 
conjugate exponent, we obtain the term 



E 



e~ kB "\^ -*p s \ kp ds 



for which we have already obtained a uniform bound, and the term 



E 



\VU n {s,c 


- W n (s, </>.) 


2s 






2s 



ds 



(11) 



For this term, we proceed as in the proof of Lemma (TH The key point is the 
estimate of the term 



E 



2TTnl „ J,n,r\|2 £ 



\V 2 U n (s, 



ds 



dr 



where 



= rx + {l-r)y+ / rb n {s, <j>») + (1 - r)b{s, s )ds + W t . 
Jo 

We can conclude as in the proof of Lemma dH if we use the result of Lemma 
IHfl In particular, (fTTj) is controlled by H^Hl 1 ?. ■ 
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We return to the proof of Lemma [3J Thanks to the uniform bounds 
obtained for the expectations in the second and third lines of fllOp . we can 
pass to the limit in n to obtain 



limsup sup E 

n ie[0,T] 



< limsup C a (\<t> {x) - Mv)\ a + f(n) a ) 



+ C a ,x lim sup E 

n 

< C a \x - y\ a + C a , x 
Using Gronwall lemma we get 



-B™ L/.n 



limsup sup E 

n te[0,s] 



limsup sup E 

n te[o,T] 



< C aX T \x - y\ 



ds. 



(12) 



We can now get rid of the exponential factor using again Holder inequality 



limsup sup E 

n te[o,T] 



\ip n - tp\ a 



< lim sup < E 

n I 

< C pXT \x - y\ a . 



e 2B T 



1/2 

sup E 

te[o,T] 



2a 



1/2 



With a = 2p, redefining as 1/2 of the process defined above and using 
the relation ([7]), we can finally transport this bound to the flows 0": 



limsup sup E 

n t€[Q,T] 



<C p limsup sup E - ip t \ p + f p (n) 

n MG[0,T] I- - 1 

< c p ,\,t\ x - y\ p ■ 



Remark that all the estimates found are uniform in x,y G B(0,R), so that 
we have obtained the desired result for the forward flows. But since the 



backward flows 



it, n i 



and 



are solution of the same SDE driven by the 



drifts — b n and —b, the same result holds for them too. ■ 
Lemma 5. For every p > 1, there exists Cd )P) T > such that 

sup sup E [|V0o' n (20| P ] < Cd, P ,T 



(13) 



uniformly in n. 
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Proof. Again, since the backward flow satisfies the same SDE of the forward 
flow with a drift of opposite sign, it is enough to show that the uniform bound 
f fl~3|) holds for the forward flows. Let 6 n and £ n be the derivatives of <j) n and if) n , 
respectively. Since 0™ = {it)' 1 °^°7o\ from fl25l) we have \8™\ p < Cd, p \$, n \ p . 
Therefore, we only need to show that the estimate f|T3l) holds for the flow ip n , 
which solves 

der(x) = wu n (t, <t> n t {x)) ${z) dt + v 2 u n (t, <k(x)) e»(x) dw t 

with initial condition £q (x) = Id. For the rest of the proof we take any fixed 
x G M. d . VU n is bounded uniformly in n and the function V 2 U is at least 
in LfQ, so that the last term is the differential of a martingale (dM t n ) due to 
Lemma [113 By Ito formula we have therefore 

d|£T < C|eri P dt+dM^+|eri p - 2 Tr {[W 2 U n (t, ff] [V 2 U n (t, #(*)) g] T )dt . 

The constant C can be chosen independently of n, and the trace of the matrix 
in the last term above can be controlled by a constant C Pt d , depending on 
p and the dimension d of the space, times |£™| 2 | V 2 £/ n (t, 0™)| 2 - Introduce the 
process 

A?:=C p , d [ \V 2 U n (s,r s )\ 2 ds . (14) 
Jo 

This is a continuous, adapted, non decreasing process, with Aq = and, 
due to Lemma [T5| E[v4J] < C uniformly in n. Lemma IT41 even provides the 
bound E[e fcA ™] < C[|i7»[| for any real constant k. We can therefore find a 
bound uniform in n reasoning as in Lemma [12j We find that 

and after integrating and taking the expected value one obtains 

Efe-^I&T] < M p + C fEie-^lCr] ds . 

Jo 

Take the supremum over all t G [0, T] and apply Gronwall inequality to get 

sup E[e- A "i&T] < ch&i p = <w , 

*G[0,T] 

uniformly in n and x G M d . Using Holder inequality as in the proof of the 
previous lemma, we finally obtain estimate f [T3|) for the derivative of the flow 
ip n , and this concludes the proof. ■ 
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3. Main Result of Existence of Weakly Differentiate Solutions 

Consider the SPDE in Stratonovich form 

du t „ „ dW n 

— + • Vu + aVu o — - = , u\ t=0 = u . 

at at 



The Ito formulation (as explained in detail also in 15j) is 

du + b ■ Vu dt + aVu dW = — Au dt , «|t=o = u o ■ 
In this section we assume b £ LL with p, g satisfying condition (OQ). 

Definition 6. Assume that b £ L^, p, g as in (Qp. VFe say iaai u is a 
weakly differentiable solution of the SPDE if 

1. w : f2 x [0, T] x M d — > R is measurable, J u (t, x) <p (x) dx (well defined 
by property 2 below) is progressively measurable for each <p £ (R d ) ; 

2. P (u (t, ■) £ n r >iW/ ' c r (R d )) = 1 /or every t £ [0, T] and both u and Vu 
are inC°({0,T];r) r > 1 L r (ttxR d )); 

3. for every <p £ Cq° and t £ [0, T], with probability one one has 



u(t, x)ip(x) dx + J J b(s, x) ■ V«(s, x)ip(x) dxds 

u(s,x) d Xi tp (x) dx ) dW l s 



= J u (x)ip(x) dx + cr^T (j 

a 2 [* f 
+ — J J u(s, x)Aif{x) dxds . 



Remark 7. The process s h- > := j u (s, x) d Xi ip (x) dx is progressively 
measurable by property 1 and satisfies J Q T |Y,?| 2 ds < oo by property 2, hence 
the ltd integral is well defined. 

Remark 8. The term J Q * J b (s, x) ■ Vn (s, x) ip (x) dxds is well defined with 
probability one because of the integrability properties in (t, x) ofb (assumptions) 
and Vu (property 2). 

Remark 9. From 3 it follows that j u (t, x) ip (x) dx has a continuous adapted 
modification, for every (p £ C£°(M. d ). 
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Let 4>t (oj) : M. d — > R d be the a-H61der continuous stochastic flow of home- 
omorphisms, for every a G (0, 1), associated to the SDE 

dXf = b (t, Xf) dt + dW t , X* = x 

constructed in [13| . The inverse of <pt will be denoted by 0q. 

Theorem 10. Assume b <E L q p withp,q as in (QJj. Ifu G n r >iW 1,r (M. d ) then 
u (t, x) := uq {4>q{x)) is a weakly differentiable solution of the SPDE. 

Proof. Step 1 (preparation). The random field (u,t,x) H- uo(<Pq(x)) is 
jointly measurable and (u>,t) (-)■ J u (<Pq(x)) (p (x) dx is progressively mea- 
surable for each ip e C^°(lR a! ). Hence part 1 of Definition O is true. We could 
prove part 2 by chain rule and Sobolev properties of </>q (x). However, a di- 
rect verification of part 3 from the formula u (t,x) := Uq (0q (x)) is difficult 
because of lack of calculus. Hence we choose to approximate u (t, x) by a 
smooth field u n (t,x); doing this, we prove both 2 and 3 by means of this 
approximation. 

Let be a sequence of smooth functions which converges to u in 
W 1,r (M> d ) and uniformly on W 1 . It is easy to check that these properties 
are satisfied for instance by Uq (x) = f 9 n (x — y) uq (y) dy when 6 n are usual 
mollifiers; for instance, the uniform convergence property comes from 

K ( x ) ~ u o (x) | < J Q n [x - y) \u (y) - u (x) I dy 

<C J 9 n (x - y) \y - x\ a dy = C J 9 n (y) \y\ a dy 

because u G C°' a . 

Let 0™ (jjj) : M d — > M. d be the stochastic flow of smooth diffeomorphisms 
associated to the equation 

dxr = b n (t, xn dt + dw t , x x ' n = x , 

where b n are smooth approximations of b as considered in the previous sec- 
tion, and let 0q" be the inverse of Then u n (t,x) := Uq (^q" (%)) is a 
smooth solution of 

du n + b n ■ V« n dt + a Vu n dW = —Au n dt , u n \ t=Q = u% , 



12 



see [18|, Theorem 6.1.5], and thus it satisfies 



u n (t, x) <f (x) dx + / / b n (s, x) ■ Vu n (s, x) (f (x) dxds 



= J UQ(x)f(x)dx + a^T^ (^J u n (s, x) d x .(f (x) dx^j dW l s 

a 2 f f f 

+ — J J u n (s, x) Af (x) dxds 

for every ip G C^°(M d ) and t G [0, T], with probability one. We need to 
establish suitable bounds on u n (t, x) and suitable convergence properties of 
u n (t, x) to u(t, x) in order to apply Lemma [[6] - which is the first step to 
obtain the regularity properties of u of point 2 of Definition [6] - and pass to 
the limit in the equation. More precisely, for every if G C^°(M. d ), t G [0, T] 
and bounded r.v. Z we have 



E 



Z u n (t, x) f (x) dx 



+ E 

d 



/ u% (x) f (x) dx + a V E Z 

J i=1 L Jo 



b n (s, x) ■ Vu n (s, x) f (x) dxds 



u n (s, x) d x xp (x) dx dW ! 



+ -E 
2 



Z J I u n (s, x) Af (x) dxds 



(15) 



We shall pass to the limit in each one of these terms. We are forced to use 
this very weak convergence due to the term 



E 



Z J J b n (s, x) ■ V« n (s, x) f (x) dxds 



where we may only use weak convergence of Vw n . 

Step 2 (convergence of u n to u). We claim that, uniformly in n and for 

every r > 1, 



sup / E 

te[or] JR d 



sup / E 

t€[0T] JRd 



u„(t,x)| ? 
Vu n (t,x)\ 



dx < C r 



dx < C r 



(16) 
(17) 
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Let us show how to prove the second bound; the first one can be obtained in 
the same way. We use the representation formula for u n and Holder inequality 
to obtain 



E 



\Vu n (t,x)\ r 



dx) < E 



\2r 



dx / E 



V0Q (x)\ 



dx . 



The last integral on the right hand side is uniformly bounded by (1131) . Also 
the other integral term can be bounded uniformly: changing variables (recall 
that all functions involved are regular) we get 



E 



|V<W B (x)) 



\2r 



dx < 



Vu n {y)\ 2r E[J m ]dy 



where J^{ y ) is the Jacobian of <fit(y)', this last term can be controlled using 
again Holder inequality, f|T3|) and the convergence of Uq in W 1,r (for every 
r > 1). Remark that all the bounds obtained are uniform in n and t. 

We consider now the problem of the convergence of u n to u. Let us first 
prove that, given t G [0, T] and <p G C ( 



oo /"itudN 



P — lim / u n (t, x) (p(x)dx = u (t, x) (p (x) dx 



(18) 



(convergence in probability). This is the first assumption of Lemma [T6l and 
allows also to pass to the limit in the first term of equation (|T5l) using the 
uniform bound (1161) and Vitali convergence theorem (we are on the compact 
support of the test function <p). Since 



u n (t, x) (p (x) dx — Uq (0o' n (x)) <P (x) dx 

(uq - u ) (0o'" {x)) V (x) dx+ I u (0o' n (x)) ip (x) dx 



using Sobolev embedding W 1,2d C 0,1//2 we have 



[u n (t, x) — u (t, x)) ip (x) dx 



< \\uq - woIIloc HvIIli 

+ / \(f>o n ( x ) -K( x ) 



1/2 



dx . 
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The first term converges to zero by the uniform convergence of Uq to u$. To 
treat the second one, recall we have proved property Hence 



lim E 



^>n n (x) - <bi (x) I dx 







and thus 



P — lim 

n—>oo 



l^o" ( x ) ~ 0o ( x ) I dx = 0. 



Property (|18|) is proved. 

Similarly, we can show that, given tp e 



P — lim 

n— >oo 



(u n (t, x) — u (t, x)) ip (x) dx 



dt = 0. 



(19) 



This implies that we can pass to the limit in the last two terms of equation 
(IT5"|) . Indeed, property f|T9l implies that 



P- lim j {^j u n (s,x) d x .<p (x) dx^j dW l s = [ ( / u (s, x) d x xp (x) dx 



10 



■J dWl 



for each i = l,...,d. Moreover, 



E 



u n (s, x) d Xi (f (x) dx J dWg 



2 










= E 


/ 


J 






Jo 





u n (s,x) d Xi ip (x) dx 



ds 



which is uniformly bounded in n due to ffTol) . By Vitali convergence theorem 
we obtain that 



lim E 

n— >oo 



Z I [ I u n (s, x) d Xi ip (x) dx ) dW l s 



= E 



u (s, x) d Xi cp (x) dx J dW ; 



The proof of convergence for the last term of equation ( JTBT) is similar. 

Step 3 (regularity of u). Let us prove property 2 of Definition [6j The 
key estimate is property (|T3"]) . 

Given r > 1 and t e [0,T], let us prove that P (u (t, ■) e W^(R d )) = 1. 
We want to use Lemma [T6l with F = u, F n = u n . Condition 1 of Lemma 
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is provided by (ITS|1 . It is clear that u n (£, •) e W 7 ^'^!^) f° r -P _a - e - w > so that 
condition 2 follows from the uniform bound on Vu n obtained in ( I17p . We 
can apply Lemma [161 and get u (£, ■) 6 H^ '^(R d ) for P-a.e. w. 

Let us prove the second part of property 2 of Definition |6j We have, from 
Lemma [T6l and (E 



E 



|Vw (£, x)\ r dx 



< limsup E 



|Vm„ (£, x)\ r dx 



Br 



< C r 



for every R > and £ G [0, T]. Hence, by monotone convergence we have 



sup E 

te[o,T] 



|Vit (£, x)| r dx 



< a 



(20) 



A similar bound can be proved for u itself: using (fT6|) . the convergence in 
probability proved in the previous step and Vitali convergence theorem we 
get that for any r' < r, R > and uniformly in time, 



E 



B r 



\u{t,x)\' 



dx = lim 



E 



Bp 



\u n (t,x)\' 



dx < C r ; 



by monotone convergence it follows that 



sup / E 

te[o,T] JR d 



\u(t, x)| r 



dx < C r 



Step 4 (passage to the limit). Finally, we have to prove that we can 
pass to the limit in equation ( !T5|) and deduce that u satisfies property 3 of 
Definition |6j We have already proved that all terms converge to the corre- 
sponding ones except for the term E Z J * J b n (s, x) ■ V«„ (s, x) (x) dxds . 

We do not want to integrate by parts, for otherwise we would have to assume 
something on div b. Since b n -> b in L q p = L q ([0,T] ; L p (R d )), it is sufficient 
to use a suitable weak convergence of Vu n to Vu. Precisely, 



E 



Z I J b n (s, x) ■ Vu n (s, x) if (x) dxds 

— EzJ J b(s,x) ■ Vn (s, x) cp (x) dxds 
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J« (t) = E 
4 2) (t) = E 



(6 n (s, x) — b (s, x) ) • Vm„ (s, x) </? (x) dxds 



Z J J if (x) b (s, x) • (Vm„ (s, x) — Vu (s, x) )dxds 



We have to prove that both J„ (i) and J n 2 ^ (t) converge to zero as n — > oo. 
By Holder inequality, 

(t) < C || 6 n — b\\ Lq ^ Q ^. LP ( Rd ^ E llVwnll^QQ^j.^p/^N 

where l/p + 1/p' = 1 and 1/g + 1/q' = 1. Thus, from (jTTJ) , In 1 -* (i) converges 
to zero. 

(2) 

Let us treat Ik (t). Using the integrability properties shown above we 
have 



E 



(p (x) b (s, x) ■ (V«„ (s, x) — Vu (s, x) ) dxds 
— J E J Zip (x) b (s, x) ■ (Vu n (s, x) — Vm (s, x) ) dx 



ds 



The function 



/in (s) := E 



y (x) b (s, x) ■ (V«„ (s, x) — Vu (s, x) ) dx 



converges to zero as n — > oo for almost every s and satisfies the assumptions 
of Vitali convergence theorem (we shall prove these two claims in Step 5 
below). Hence 1^ (i) converges to zero. 

Now we may pass to the limit in equation (|T5|) and get 



E 



Z u(t,x) ip (x) dx 



+ E 
d 



b (s, x) • Vw (s, x) ip (x) dxds 



y Mq ( x ) V 9 ( x ) dx + a E, Z J yJu(s,x)d Xt <p(x) dxj dW t 



+ -E 
2 



Z ^ / m (s, x) Ay> (x) dxds 



The arbitrariness of Z implies property 3 of Definition El 
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Step 5 (auxiliary facts). We have to prove the two properties of h n (s) 
claimed in Step 4. Recall we may use Lemma [TH] at each value of time. It 
gives us 



E 



d Xi u (s, x) (p (x) Z dx 



lim E 

n— >oo 



d Xi u n (s,x) ip (x)Zdx 



(21) 



for every ip G C^°(M. d ) and bounded r.v. Z, at each s G [0,T]. We have 
beL q ([0,7] ;U>(R d )), hence b (s, •) G L p (R d ) for a.e. s G [0,7). The space 
Co°(R d ) is dense in L p (M. d ). We may extend the convergence property (j2ip 
to all G L p (R d ) by means of the bounds (JTZD and (j20|) . Hence fc, n (s) -»■ 
as n — )■ oo, for a.e. s G [0, 7]. 

Moreover, for every e > there is a constant Cz,<p,e such that 



/ ^ +e (s)ds < C^, e [e\ [\b(s,x)\ 1+£ - ( | V«„(s, x) 1 1+£ + \Vu(s,x)\ 1+£ ) dx 
Jo Jo L Vbj? v ' 



<C^, £ ||6|m e |E 



ds 



l+e 



Vw n is, x) dxds 



+ C^, £ ||6|g £ |E 



^B fl 
T 



Vw (s, x) dxds 



^B fl 



for a suitable r depending on e (we have used Holder inequality). The bounds 
(|T7|) and (I20I) imply that h n +e (s) ds is uniformly bounded. Hence Vitali 

theorem can be applied to prove that I ( n ] it) = $1 h n (s) ds -> as n — > oo. 
The proof is complete. ■ 



4. Uniqueness of Weakly Differentiable Solutions 
Theorem 11. Weak solutions of Definition^ are unique. 

Proof. Let u % be two weakly differentiable solutions of equation 



du , „ „ dW n 
— + b ■ Vu + ctVm o — — = 
at dt 



U\ t =o = U 



Then u := u 1 — u 2 is a weakly differentiable solution of 

u\ t =o = 



— + 6 ■ Vu + ctVm o — — = . 
at dt 



(22) 
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We want to prove that u is identically zero. We divide the proof in three 
steps. 

Step 1 (Equation for u 2 ) The first step consists in proving that u 2 is also 
a weakly differentiable solution of 



du 2 , „ 2 „ 2 dW 
— + b- Vu 2 + aVu 2 o— — = 0, u\ t= o = 
at at 

namely that 

J u 2 (t, x) if (x) dx + J J b (s, x) ■ Vu 2 (s, x) <p (x) dxds 

= ° E jf (/ u2 ( s ' x ) w dx ) dw * 



(23) 



i=i 

2 /•* 



+ — y y -u 2 (s, rr) Ay? (x) dxds 

for any </? G C^°(lR d ). Let # e be a sequence of standard mollifiers. From the 
definition of weak solution, using y> e y {x) = 9 £ (y — x), we have 

Ue(t,y)+ b(s,y) -Vu £ (s,y)ds 
Jo 

+ f\v(s,y)odWi= [ t R £ (s,y)ds, 

i=1 Jo Jo 

R £ (s,y)= J (b(s,y) - b(s,x)^Vu(s,x)9 £ (x - y)dx 

The function u £ is smooth in space. For any fixed y, by Ito formula we have 
du 2 (t,y) = 2u E (t,y) du £ (t,y) 

d 

= - 2u £ (t, y) b(t, y)Vu e (t, y) dt - 2a u £ (t, y) ^ d yi u £ {s, y) o dW l s 
+ 2u £ (t,y)R £ (t,y) dt 



i=i 
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which, rewritten in the week formulation using a generic test function ip, 
reads 



u 2 £ (t,y)tp(y)dy + J J b(s,y)Vu 2 (s,y)tp(y) dyds 

+ a J2f{f yMv) d ^) ° dw * = If 2u ^ s ' y) R ^ s ' vMv) d y ds ■ 

We want now to pass to the limit for e — >• in the different terms. Since for 
every t, u £ — > u uniformly on compact sets, by dominated convergence the 
first term tends to 

u 2 (t,y)ip(y) dy. 



For the following terms, we consider s fixed. Using Holder inequality and 
the convergence of ||Vw £ ||lp — > ||Vw||lp on compact sets (recall that ip is of 
compact support) for any p > 1, we have 



f b(s, y)<p(v) (V«e(s, y) - Vm 2 (s, y)^dy 



< || & (s,2/)^(2/)|| ir /||V^(s,?/) - Vu 2 (s,y)\\ Lr 

< C r Ml p \\Vu 2 e (s,y) - Vu 2 (s,y)\\ r , , — >• , 



which is enough to obtain the convergence of the second term. In the same 
way one obtains also the convergence of the third term. As for the term con- 
taining the commutator R e , we can use again Holder inequality, the uniform 
convergence of u £ , the equi-boundedness of Vu £ in LP for every p > 1, and 
the continuity in mean (for a.e. y) of the function b G L p (M. d ). This proves 

(EHJ). 

Step 2 (equation for v 2 ) We have that u is a.s. continuous in space and 
time (and therefore locally bounded) and by definition of weak solution Vw G 
L r ([0, T] x R d ) for every r > 1 a.s.. It follows that f(s) = J Vm 2 (s, y)<p(y) dy 
is still a.s. a function in L r (0,T). This means that, writing ( )23|) in Ito form, 
the stochastic integral is a martingale and 

E[-u 2 (t,x)] ip(x) dx + J J b(s, x) ■ VE [u 2 (s, x)] (f(x) dxds 
Y f f^[u 2 (s, x)] Aip(x) dxds. 
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Hence v(t, x) = E [u 2 (t, x)] satisfies 



v(t, x)f{x) dx + / / 6(s, x) ■ Vv(t, x)(p(x) dxds 



a 2 ft f 

= — J J v(t, x)A(p(x) dxds 

and is fairly regular: v G C°([0, T]; W 1 ' r (W 1 )) for r > 1. This follows by 
Holder inequality because 

/ \Vv(t,x)\ r dx = j \K[uVu]\ r dx < J (E[M 2 ]E[|Vu| 2 ]y /2 da; 

< ( J E[M 2r ] V) 1/2 ( J E[|Vn| 2r ] 2 dx) 1/2 < C , 



uniformly in t (similar computations provide the same result for the function 

Thanks to its global integrability properties, using approximating func- 
tions as in the first step, one can prove that v solves 



/ 



v 2 (t, x) dx + (7 Z 



-2 



|Vf (t, x)\ 2 dxds 



b (s, x) ■ Vf (t, x) v (t, x) dxds . (24) 



Step 3 (final estimates) We want to find suitable bounds on the last term 
of fl24l) allowing to apply Gronwall inequality. This will complete the proof. 
For every t G [0, T], we have 



v b ■ Vf da; 



< 



|Vf| dx 



1/2 



v 2 16| 2 dx 



1/2 



J v 2 \b\ 2 dx< (J v 2r dx^ ' (j \b\ p dx^j 



\ 2/p 



where 1/r + 2/p — 1 namely 1/r = 1 — 2/p = (p — 2)/p : 

p 

r = . 

p- 2 
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One has the interpolation inequality 



l/ot / f \ 1—s / f \ s .-, 



i/Mx ) < ( / «Mi ] [ I \Vv\ z dx) , S = ^2a~ d ' 



The idea of the result comes from: W s ' C L a for 1/a = 1/2 — s/d, 
namely 4 = \ — - — ^p^, s = ^^d; and then 

l/a 



v a dx } < \\v\\ Ws ,2 < ||w||^2 S ||W||^/1,2 



Let us put everything together: 



l/r / p x 2/a / n \ 1 ~ S f f \ S 

v 2r dx) = / v a dx) < / v 2 dx) / \Vv\ 2 dx 



r = — - — , a = 2r , s = — d 

p - 2 ' ' 2a 

namely 

_2r-2 ,2^-2 d 

S Caj a 7} 

p— 2 r 

Thus we have proved: 

J v 2 \b\ 2 dx< (j v 2 dx'j " (^J \Vv\ 2 d x y (^J \b\ p dx^j ' 
and we can bound the last term in 



v b ■ Vf dx 



< / |Vu| dar 



p-\-d p — d 

f dx 



\b\ p dx 



i/p 



Recall that ab < ^ + £ ± + 1 = 1. Then, with s = r = we have 



2p 



2p 



v b • Vf dx 
Therefore 



< — / | Vvr dx + C / Wx / |6| p dx 



2 

p—d 



v 2 (t,x)dx<C [ v 2 dx / |6| p dx ds 



2 

p—d 
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hence we may apply Gronwall lemma and deduce j v 2 (t, x) dx = if 



T ' 



p — d 



|6| p da; ds < oo. 



We know that 



T 

\b\ p dx) ds < oo 



for certain p, q > 2 such that - + - < 1 . Then 

2 <q 



p — d p 

because 2p < qp — qd, - < 1 — -. The proof is complete. ■ 

5. Appendix: Technical Lemmas 

For completeness, we collect here some modifications of known results 
used in Section |2j We will use the notation introduced there. 

Lemma 12. Let U n be the solution of the PDE for f = b = b n , as defined 
in Section^ Then 

i) U n (t,x) and VU n (t,x) converge pointwise in (t,x) to U(t,x) and 
W (t, x) respectively, and the convergence is uniform on compact sets; 

ii) there exists a A for which sup ta; \VU n (t,x)\ < 1/2; 
Hi) \\V 2 U n (t,x)\\ LUT) <C . 



Proof. The result of the second point is proved in 13j, Lemma 3.4] for a 
fixed n, but inspecting the proof we see that all the bounds obtained depend 
on || & || £«, but never on b itself. Since ||& n ||L« — * the uniformity in n 

follows. 

To prove the other two points, set V n := U n — U; then 

d t V n + ^AV n + b ■ W n - XV n = -{b n -b)- (Id + W n ) , V n (T, x) = . 

From the bound ([5]) on the solution provided by Theorem El we obtain 

IMI^ < iV||&^ - ftllxi ^ . 



'2,p 
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It follows that U n — > U in #f jP , which proves the last point. Since by [17 . 
Lemma 10.2] U, U n , VU and VU n are all Holder continuous functions, there 
exists a subsequence (that we still call U n ) s.t. U n — > U and VU n — > VU for 
every (t, x) and uniformly on compacts. ■ 

Remark 13. The following results hold uniformly in n because, as remarked 
in the previous proof, all the bounds obtained depend on the norm ofb. 

i) From fl3 . Lemma 3.5] we have 

sup sup <2; (25) 

n te[0,T] v ; 

ii) from the uniform boundedness of the coefficients (U n and VU n ) of the 
SDE we get 



sup E 

te[o,T] 



i/j?(x)\ a < C(l + |x| a ) . (26) 



Lemma 14. For every n, both the process A n defined by and the one de- 
fined by ( fl^y are continuous, adapted, nondecreasing, with Aq = 0, K[Aj] < 
C and for every feel, E[e kA T] < C. The constant C can be chosen inde- 
pendently of n. 

Proof. For the process defined by (Q the proof follows the same steps of the 
proof of 0, Lemma 7]. We only remark that the function U is the solution 
of a different PDE, but it has the same properties in terms of regularity. 
Moreover, the flows and <p n solve two SDEs with different drifts b and b n , 



which means that in the proof one has to use twice the result of [12l . Corollary 
13], once for every drift. 



For the process defined by (fT4|) . the result is already contained in [12, 
Corollary 13]. ■ 

Lemma 15. Let f n be a sequence of vector fields belonging to L^, convergent 
to f E L^. Then, there exists e > 1 s.t. 

E\ f \f n (s,^)\ 2£ ds] < C < oo . (27) 

Proof. To prove the result for a fixed n one can use flil . Corollary 13] and 
follow the proof of [l2j, Lemma 8 and Corollary 9], which still works due to 
the strict inequality in the conditions imposed on p, q. Then, since all the 
bounds only depend on the norm of / but never on the function itself, one 
obtains that ( 1271) is uniform in n. ■ 
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6. Appendix: Sobolev Regularity of Random Fields 

Let r > 1 be given. We recall that / G (R d ) if / G U loc (R d ) and 
there exist gi G U loc i = 1, such that 



/ (x) d x .(p (x) dx = - / ft (x) (a;) cfe 



for all if G (]R d ). When this happens, we set d x J (x) = gi (x). From the 
definition and easy arguments one has the following criterion: if / G L\ oc (lR d ) 
and there exist a sequence {/„} C (R d ) such that f n — > f in L\ oc (lR d ) 
(or even in distributions) and for all R > one has a constant (7r > 
such that / |V/„ {x)\ r dx < C R uniformly in n, then / G W,^ r (R d ). This 
criterion will not be used below; it is only stated for comparison with the 
next result. 

Let now F : Q x R d — > R be a random field. When we use below this 
name we always assume it is jointly measurable. 

Lemma 16. Assume that F (u, •) G U loc (R d ) for P -a.e. lu and there exist 
a sequence {F n } neN of random fields such that 

1. F n (uu, •) — > F (uu, •) in distributions in probability, namely 

P — lim / F n (u>, x) if (x) dx — / F (u, x) (f (x) dx 

n ^°°jRd J R d 

for every ip G (R d ) ; 

2. F n (u, •) G Wfal (R d ) for P-a.e. u and for every R > there exists a 



constant Cr > such that 



E 



uniformly in n. 
ThenF(u; r )eWl; r c 



\VF n (x)\ r dx 



for P-a.e. u, 



< C 



R 



E 



/ d Xi F (•, x) (f (x) Z dx = lim E / d Xi F n (•, x) p (x) Z dx 



(28) 
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lim 

n— >oo 



for all (p G Cq 
P 

for all ip G Co 
E 



and bounded r.v. Z , 
d Xi F n (u,x) <p (x) dx = - 
, and for every R > 



F(u,x)d Xi ip(x)dx (29) 



IVF (x)\ r dx 



< lim sup E 

n— >oo 



\VF n (x)\ r dx 



(30) 



Proof. Given R > 0, there is a subsequence {n k } and a vector valued 
random field G such that VF„ fc converges weakly to G in L r (Q x as 
— ?> oo. Taking i? G N, we may apply a diagonal procedure and find a single 
subsequence {n^.} and vector valued random fields G R , R G N, such that 
VF nk converges weakly to G R in U (Q x as k — )■ oo, for each i? G N. 
Using suitable test functions, one can see that G R ' = G R on SI x B R if R' > R. 
Hence we have found a single vector valued random field G, such that V_F ni , 
converges weakly to G in U (Q x Br), as /c — > oo, for each R G N and thus 
for each real i? > 0. At the end of the proof, G will be identified by VF, 
independently of the subsequence {n^}. Thus, a fortiori, the full sequence 
VF n converges weakly to G in U (fi x R d ), as n — > oo. For this reason, to 
simplify notations, we omit the notation of the subsequence. 
For each (p G C 3O (]R d ), by assumptions 1 and 2 



/ 



F (u>, x) d Xi (p (x) dx — lim / F n {u, x) d Xi tp (x) dx 

= - lim / d Xi F n (u,x)(p(x)dx , 



the limits being understood in probability. For each bounded r.v. Z, this 
implies that 

lim / d Xi F n (u, x) (p (x) Z (cu) dx = — / F (u>, x) d Xi (p (x) Z (u>) dx 

n->ooJ Rd ' JRd 

in probability. This limit holds also in L 1 (Q) by Vitali convergence criterion 
because, by Holder inequality, 

E / d Xi F n (to, x) if (x) Z (u) dx 



< Cr, p II^IIloo ii^iu°°e 



/ \d Xi F n (u,x)\ p dx 

J B R 
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uniformly in n, for some p > 1, and with R such that Br contains the support 
of ip. 

From the weak convergence above, we also get that 



lim E 

n— >oo 



/ d Xi F n (-,x)tp(x)Zdx 



E 



Gi (•, x) tp (x) Zdx 



Hence 



E 



Gi (■,x)<p (x) Zdx 



-E 



[ F(-,x)d Xi ip(x)Zdx 
JR d 



By the arbitrariness of Z this gives us 



F(u,x) d Xi ip (x) dx 



Gi (cj, x) <p (x) dx 



(31) 



for P-a.e. u. This is the identification of G mentioned above, which implies 
the weak convergence of the full sequence VP n . 

Identity (l3Tj) holds P-a.s. for every ip G C^°(M. d ) a priori given. Thus it 
holds P-a.s., uniformly on a dense countable set V of test functions ip, dense 
for instance in W l( £ (IR d ), 7 + ^7 — 1- Using the integrability properties of 
F (u, •) and Gi (u, •) we may extend identity §T§ to all <p G (M d ) and 
thus all <p G C^°(lR d ), uniformly with respect to the good set of u for which 
it holds on T>. 

Thus, from identity fl3Tj) in the stronger form just explained, we deduce 



by definition - that F(u,-) G Wj£(R d ) for P-a.e. u. And VP (u, x) 



G (u, x). We immediately have (gg) and (j2gj) . 

We have shown that, for every function £(u,x) of the form 

m 

£(u,x) = ^(f k (x)Z k (u) , 



fc=l 



with ip and Z as above, we have 



E 



VP(-,x)£(-,x) dx 



Br 



< lim 

n— >oo 



E 



VP n (-,x)£(-,x) dx 



< E 



|£(-,x)| r dx 



l/r' 



lim sup E 



ivp„(-,x)rdx 



Ub r 



l/r 
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(in the last passage we have used Holder inequality). The set of functions £ 
introduced is dense in U (fl x B R ), so that 



E 



Br 



|VF(-,a;)| r da; 



l/r 



VF 



< sup 

U\\ L r><l 



< lim sup E 



L r (QxB R ) 

E 



/ VF{;x)£(;x)dx 

J Br 



[ \VF n (.,x)\ r dx 

J Br 



l/r 



This completes the proof. 
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